Spin dynamics of stripes 
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The spin dynamics of stripes in high-temperature superconductors and related compounds is 
studied in the framework of a spin-wave theory for a simple spin-only model. The magnon dispersion 
relation and the magnetic structure factor are calculated for diagonal and vertical stripes. Acoustical 
as well as optical bands are included in the analysis. The incommensurability and the ir resonance 
appear as complementary features of the band structure at different energy scales. The dependence 
of spin- wave velocities and resonance frequencies on the stripe spacing and coupling is calculated. 
At low doping, the resonance frequency is found to scale roughly inversely proportional to the stripe 
spacing. The favorable comparison of the results with experimental data suggests that the spin-only 
model provides a suitable and simple basis for calculating and understanding the spin dynamics of 
stripes. 

PACS numbers: 75.10.Jm, 74.72.-h, 75.30.Fv, 76.50.+g 



I. INTRODUCTION 



The evidence for the formation of stripes in high- 
temperature superconductors (HTSC) and related ma- 
terials increases continuously. After the theoretical 
predictionErBu of stripes as a combined charge and 
spin-density wave phenomenon, years passed until a 
broad interest was triggered by experiments on insu- 
lating La2_a;Sr x Ni04 +( 5 1LSNO) and superconducting 
La2- x Sr-CjiiQ4 (LSCO).eEI More recent experimental 
evidenceoQtm for stripes in the paradigmatic HTSCs 
YBaCu0 6+ 5 (YBCO) and Bi 2 Sr 2 CaCu 2 8+( 5 (BSCCO) 
strengthens the expectation that stripe formation in 
doped layered perovskites is quite generic. 

In spite of the striking evidence for stripes in HTSCs, 
the causal connection between stripe formation and su- 
perconductivity still is a mystery. It is puzzling that both 
phenomena coexist and that, neyertheless, stripes tend to 
suppress superconductivity.E30 For this interplay spin 
order is more relevant than charge order. In particular, 
the strength of spin fluctuations appears to play a central 
role. Static spin order seem to be much less compatible 
with superconductivity than dynamic spin order. 

At present, one important open question is to what 
extent the stripe model can account for spin fluctua- 
tions not only at low energies, where the incommensu- 
rate response is observed, but also over a wider energy 
range, including the resonance phenomenon at the anti- 
ferromagnetic wave vector (see Refs. |l^,|l3| and references 
therein) . The specific form of the dynamic magnetic re- 
sponse - including incommensurability and 7r resonance 
- gave rise tcuioubts that it could be consistent with the. 
stripe modeLfca On the other hand, there are proposalsEJ 
that both features may be rooted in a stripe-like spin- 
density wave. 

In this papst,we complement the spin- wave analysis by 
Batista et al.J— I There, the emphasis was put on the in- 
commensurate ratio between the spin spacing a and the 
stripe spacing d which gives rise to a continuous excita- 
tion spectrum. However, in many cases of interest, this 



ratio p := d/a is very close to a rational or even inte- 
ger value. Within the stripe model one actually expects 
that integer values are very stable due to a lock-in of the 
superstructure into the atomic structure. This pinning 
mechanism is considered as the origin of the so called 
T/8 conundrum' in the cupratesfJ i.e., the stability of 
p = 4 over a considerable doping range. Detailed mea- 
surements of the spin-excitation spectrum are available, 
close to integer p: p =f-3 in LSNOP p = 4 in LSCO,Ell 
and p = 4E3 and p = 3ii in YBCO. 

In order to test whether these experiments can be con- 
sistent with the spin-wave excitation spectrum of a stripe 
model in the simplest and most transparent case, we 
therefore examine integer p. Particular attention is paid 
to the spin-wave band structure in the vicinity of the 
antifcrromagnetic wave vector. While the incommensu- 
rability as zero-frequency response is fixed by the geome- 
try of the model, we calculate the spin- wave velocities at 
the incommensurability and the tt resonance as dynamic 
features. We evaluate the dependence of these quanti- 
ties on the incommensurability (respectively, the doping 
level) and the exchange coupling across the stripes. By a 
quantitative comparison, we determine the value of the 
exchange coupling across the stripes as the only a pri- 
ory unknown model parameter. In particular, the de- 
pendence of the 7r resonance on doping is found to be 
consistent with experiments. 

Our course starts in Section || with the introduction 
of the spin-only model that constitutes the basis of our 



study. The linear spin- wave theory is outlined in Sec. Ill 



In Sec. IV we present numerical results for the magnon 



dispersion relation, spin-wave velocities, tt resonance and 
the structure factor. In Sec. [V| the results of our theory 
are discussed and compared to experimental data. 



II. MODEL 

In the cuprates as well as in the nickelates, the metallic 
spins are located on square lattices in weakly coupled lay- 
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ers. Since the the interlayer coupling generally is much 
smaller than the intralayer coupling, we focus on a sin- 
gle layer. For simplicity, the holes induced by doping 
are assumed to form site-centered rivers that act like aru 
tiphase boundaries for the antiferromagnetic domains.u 
The rivers are assumed to be only one lattice spacing 
wide (cf. Fig. |). 
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FIG. 1: Illustration of vertical and diagonal stripe patterns 
with spacings p — 3 and p = 4. The hole positions are indi- 
cated by open circles and the electron positions by grey ones. 
The arrows correspond to the spin orientations in the clas- 
sical ground state. Parallelogramms outline magnetic unit 
cells spanned by A' 1 - 1 and A' 2 '. In our model, we assume an- 
tiferromagnetic exchange couplings of strength J within the 
domains (dashed lines) and A J across stripes (zig-zag lines). 

Since stripes are vertical in cuprates for doping con- 
centrations where superconductivity occurs and diagonal 
in nickelates, we study both orientations with arbitrary 
integer stripe spacing p. Furthermore, since charge order 
seems to.be static up to high temperatures, in YBCO up 
to 300K,B holes can be considered as immobile at low 
temperatures. Our analysis is restricted to T = 0. 

We are interested in collective excitations around a 
ground state, which - for classical spins - could be rep- 
resented by S(r) = 5{0,0,tr(r)} with a = ±1 on the 
electron positions and a = on the hole positions (as 
illustrated in Fig. [I]). Denoting by and A^ 2 ^ the 

primitive basis vectors of the magnetic unit cell and by 
A = miA^ 1 ' +m2A' 2 ' an arbitrary magnetic lattice vec- 
tor, the classical spin variables obey the translational 
symmetry tr(r) = cr(r + A). By placing the origin at 
a hole position we obtain the additional reflection sym- 
metry <j(r) = — <r(— r). 

For a paradigmatic and minimalistic description of 
magnetic quantum fluctuations we use a spin-only model 
with pair exchange. More cpmelicated exchange pro- 
cesses such as cyclic exchanget£lcJ may be important for 
quantitative purposes but are ignored here for simplic- 



ity. We use a generalized Jieisenberg model on the two- 
dimensional square latticed 



n 



-]T./(r,r')S(r)S(r'), 



(1) 



where the primed sums run over all spin positions with 
a 7^ 0. The exchange couplings obey the symmetry rela- 
tions 



J(r,r') = J(r',r), 

J(r,r') = J(r + A,r'+A), 

J(r,r') = J(-r,-r'). 



(2a) 
(2b) 
(2c) 



In fact, the exchange couplings may have a higher sym- 
metry corresponding to the hole lattice, which however 
will not be needed explicitly in the further analysis. To 
implement that the hole strings act as antiphase bound- 
aries between antiferromagnetic domains, we assume that 
J(r, r') = J > for nearest neighbors r, r' within the do- 
mains and J(r,r') = A J > for nearest neighbors across 
a string. 

While it is natural to assume that J should be compa- 
rable to the exchange coupling in the undoped material, 
the coupling AJ may deviate significantly. To keep the 
number of parameters small, we ignore that the exchange 
coupling even within an antiferromagnetic domain should 
depend on the position of the pair relative to the hole 
strings. 



III. SPIN- WAVE THEORY 

We address the spin dynamics in the framework of 
linear spin-wave theory (for a review in the context of 
cuprates, see e.g. Ref. 



22). In the following analytic 



part we keep the general form of the model and special- 
ize to specific stripe configurations later in Sec. |V| when 
we numerically evaluate the results of this section. From 
now on all lengths are expressed in units of the square- 
lattice spacing a. 



A. Holstein-Primakoff representation 



In the first step we flip all spins on one sublattice by 



S x (r) = a 2 {r)S x (r), 
S y (r) = a(r)S y (r), 
S z (r) = a(r)S z (r). 



(3a) 
(3b) 
(3c) 



This transformation preserves the spin commutator rela- 
tions. Thereby, we allow S to have spin S also at the hole 
sites. Although this introduces certain modes of zero en- 
ergy, as we will discuss below, it is advantageous to use 
a S with a homogeneous ferromagnetic ground state. 
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The corresponding transformed Hamiltonian reads 
H = ^J(r,r')[5 z (r)^(r') + ^(r)^(r') 

r,r ; 

+a{v)a{v')S x (v)S x {v')\ , (4) 



where we have defined the new couplings J(r,r') := 
J(r, r'W(r)cr(r') which obey the same symmetry rela- 
tions (g) as J. 

In the next step we represent the spin operators by the 
usual Holstein-Primakoff (HP) bosons via 



S + = V2S -hb, 
S~ = tfy/2S-h, 
S z = -n + S, 



(5a) 
(5b) 
(5c) 



with S = S x ± iS v . The eigenstates of the number 
operator n = b^b are restricted to n < 2S and the 
HP-operators fulfill the canonical commutator relations 
[6,6^] = 1. The linearized spin- wave Hamiltonian H sw is 
given by the terms quadratic in the bosonic operators, 



J2{fM [&t(r)6(r')+6(r)6t(r')] 



r,r' 



?(r,r') [b(r)b(r') + b*(rp(r')]}, (6s 



/(r,r') = lj(r,r')Kr)a(r') + l], 



g(r,r') = 1 J(r, r') [a(r)a(r') - 1] 



(6b) 
(6c) 



Obviously the functions / and g again satisfy the sym- 
metry relations (||). 

For further manipulations it is useful to decompose a 
vector r = A + a on the square lattice into a vector 
A = toiA^ 1 ) + m2A' 2 ' on the magnetic lattice and a 
decoration vector a. The number of vectors a is denoted 
by N (the area of the magnetic unit cell). In momen- 
tum space, the reciprocal magnetic basis Q*- 1 -* defines the 
corresponding magnetic Brillouin zone (BZ). Wave vec- 
tors k can be uniquely decomposed into k = Q + q with 
q e BZ and Q = miQ^ +m 2 Q (2) . Within the Brillouin 
zone of the square lattice there are N vectors Q which 
we denote by Q, y . 

We Fourier transform the bosonic operators via 6(r) = 
f k exp(ikr)&(k), where J k = (27r)~ 2 f d 2 k and the k in- 
tegrals run over the Brillouin zone of the square lattice 
with an area (2tt) 2 . Using these decompositions and the 
Poisson sum formula 



E^ kA - ^E*( k +Q) 

A Q 



(7) 



we rewrite the spin-wave Hamiltonian as 

+&-q-Q„& t _ q _ Q J 



-&-q-Q„&q+Q„ 



(8) 



where 



JV(q) = -y^^/(a + A,a') 

A a, a' 

x cos [qA + q(a - a') + Q,,a - Q„a'] (9) 



is essentially the Fourier transform of /, 

jjf(Qu + q, Q„ + q') = <5(q + q')^'(q)- (10) 
Analogous expressions relate G to g. 



B. Bogoliubov Transformation 

To diagonalize the Hamiltonian, we express the bosonic 
operators by canonical coordinate and momentum oper- 
ators $„(q) := $(q + Q„) and n„(q) := n(q + Q„) via 
the relations 



n,(q) = (&_ q -Q„ - ft* 



q+Q t 



(11a) 
(lib) 



In terms of these operators, the spin-wave Hamiltonian 
reads 

H ™ = \ J E{ n ^) M ~V(q)iMq) 

+^M)K^,(q)^,(q)\, (12) 



with the inverse mass matrix M _1 = F — G and the 
coupling matrix K = F + G. As a result of the invari- 
ance of the Hamiltonian under the replacement S x (r) — > 
a(r)S x (r), S y (r) — > a(r)S v (r) one can easily derive the 
symmetry conditions 



K = o-M V, 
M" 1 = aKa. 



(13a) 
(13b) 



where we have introduced the hermitian matrix a v y := 
^ a e~ 4(Ql,_Q "' )a o'(a). To simplify notation, we sup- 
press arguments q which may be considered as fixed dur- 
ing the diagonalization in v space and use the pseudo 
Dirac notation |4) := |n)) := EM") 
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with the Cartesian basis v — 1,. . . ,N. After per- 
forming the canonical transformation |<£>)) = M^ 1 / 2 ! <!>)), 
Ill)) = M^ 2 \U)) the Hamiltonian can be rewritten as 



: i J |((fi|n)) + ((ilM-^KM- 1 / 2 ^))!, (14) 

and we still have to diagonalize M" 1 / 2 KM^ 1 / 2 = W 2 
with hermitian W := M _1 / 2 <tM -1 / 2 . Introducing an 
orthonormal eigenbasis {|a)), a = 1, . . . N} of this matrix, 
W|a)) = £,a\cx}, and defining w a := |£ a | we can transform 
to normal coordinates 



-1/2, 



n„ 



E-« /2 < 



Qf|i/))n , 



and obtain 



(15a) 
(15b) 

(16) 



Transforming back to corresponding bosonic operators 
$„(q) = ^[6„(q) + 6L(-q)l. &a(q) = ^[6a(-q) - 
(l)] we obtain the final diagonal bosonic representation 
of the spin-wave Hamiltonian 



n sw = Y.J ^(q){^ + ^(q)Uq) 



(17) 



Thus, as the result of the above diagonalization we obtain 
w Q (q) as the magnon dispersion relation with the band 
index a. 

We would like to remark that the \vj) space contains 
a common subspace of eigenvectors of the matrices cr, 
M _1 and K with vanishing eigenvalues. This subspace 
is h dimensional, where h is the number of holes in the 
magnetic unit cell. These zero modes are an artifact of 
the introduction of spins S on the hole sites. All above 
manipulations, including e.g. the calculation of M 1 / 2 
and w" 1 , are welldefined on the orthogonal subspace of 
physical spins. 



C. Structure Factor 

In this section we proceed to calculate the zero- 
temperature structure factor 

5(k, W ) := E E l<F|^(k)|0)| 2 ^-c F ).(18) 

Here, |0) denotes the ground state (magnon vacuum) 
characterized by 6 Q (q)|0) = and we consider only 
single-magnon states |F) with excitation energy wp := 
E F — E . Since 

S z (k) = sX*(it-Q,'MQ^) 

v' 

-E CT (Q"') / tf(k")b(k-Q u ,+k") (19) 



with cr(Q) := J2 a e~ lCJa a'(a) preserves the magnon 
number, it contributes only to the elastic part of the 
structure factor, 

S c \k,uj) cx S ,2 ^ ( 5(k-Q)|cr(Q)| 2 (5(w). (20) 
Q 

To calculate the inelastic part of the structure factor 
(which has contributions of order S only from j = x, y) 
we express these spin components by the bosonic oper- 
ators using the transformations derived in section III B| . 



s x (q + Q„) « \/5E^ (2) (Q--Q-')^'(q) 

v' 

= /|e^ (2) (q,-q,'K- 1/2 

x((^|M- 1 / 2 | a ))[6 Q (q)+&t Q (_ q) ] 5 (2l a) 
S y (q + Q v ) « V5^a(Q,-Q^)n^(q) 

v' 

= ^E^-Q^vy 2 

a.v' 

x(^\M^\a))[b a ( q )-bi(^)}, (21b) 

where we have defined := jf^2 gL e~ lCla a 2 (a). 

Since the contributing final states are just given by the 
one-magnon states |F) = &a(q)|0) it is easy to calculate 
the inelastic part of the structure factor. Using the re- 
lations <T 2 M-^ 2 \a) =o-M 1 / 2 W|a) = ^erM 1 / 2 ^)) and 
cr 2 M-^ 2 \a} = M~^ 2 \a} we obtain 

S in (q+Q„^) = S^5 a (q + Q I ,)<5(w-Wa(q)),(22a) 

a 

S„(q + Q„) = (vlM- 1 / 2 ^)) — (alM- 1 / 2 ^)). (22b) 

At this point it may be helpful to remind that q is an 
implicit argument of u) a , M" 1 / 2 and |a)). The periodicity 
w Q (q) = w Q (q + Q) of the eigenfrequencies is absent in 
the structure factor since the coupling of an external field 
to a spin wave wave vector k = q + Q depends on Q. 

IV. RESULTS 

We now evaluate the above general analytic expres- 
sions for the magnon dispersion and the structure factor. 
Thereby we focus on our minimalistic model (cf. Sec. |J) 
with stripe spacings p = 3, 4, and 5, since these values 
correspond to doping concentrations in various experi- 
mental works as mentioned in the introduction. The ex- 
plicit comparison to experiments is postponed to Sec. [v|. 

For later reference, we briefly recall that for the un- 
doped two-dimensional antiferromagnet (which is recov- 
ered by our model in the limit p — > oo) the spin- wave 
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dispersion is given by 



WAF(k) = 2JSV 4 - [cos(2ttF) + cos(2nK)Y 



(23) 



[From now on, we refer to wave vectors k = (if, K) in 
units of 2n/a\. It vanishes at the antiferromagnetic wave 
vector kAF = {h, 5), where the structure factor shows 
maximal intensity To leading order in <5q = k — k^F) 
the low energy spin- wave excitations are characterized by 
an isotropic dispersion waf ~ wafI^qI with a spin- wave 
velocity vaf = V&JSa. 



A. Vertical case 

For vertical stripes a possible magnetic unit cell is 
given by the basis vectors = (0, 2) and A^ 2 ) = (p, 0) 
for odd or A*- 2 ) = (p, 1) for even p. Therefore we have 
N = 2p lattice sites per unit cell (cf. Fig. [l]) and 2p eigen- 
values w Q (q). Two of them (corresponding to the number 
of holes) vanish identically and we obtain p — 1 twofold 
degenerate physical bands. This degeneracy results from 
the equivalence of the two sublattices. 

The lowest, acoustical band has zeros at the magnetic 
superstructure lattice vectors. Within the Brillouin zone 
of the square lattice (we choose < H,K < 1), the 
vectors are located at (^,0) and |) for odd p or 

(p + 2p' I) f° r even P ( w ith < j < p). In the upper 
row of Fig. ||, we show u>(k) for the acoustical band as a 
density plot, where black corresponds to lu — and white 
to the upper band edge. 



p = 3 



p=4 



p = 5 



» • • ( 
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FIG. 2: Acoustical band for vertical stripes with spacings 
p = 3, 4, 5 and A = 0.5. The upper row shows density plots 
of the spin-wave dispersion, where dark regions correspond to 
low energy values. The lower row shows the acoustical band 
in the (H, K, oj) space including the weight of the inelastic 
structure factor, where larger weight corresponds to darker 
points with larger size. 

Although the dispersion relation obeys the symmetry 
Lu(k) = u(k+Q) corresponding to the period of the mag- 



netic unit cell, this symmetry is absent in the structure 
factor. In the lower row of Fig. || the acoustic band is 
rcplotted in the (H, K, lo) space using darker and thicker 
dots for points with larger values of the structure fac- 
tor (|2^). In agreement with experiments, the weights 
are concentrated near the lowest harmonic "incommen- 
surate" wave vectors Q — {\ ± ^). Higher harmonics 
of the superstructure are much weaker as already noticed 
in Ref. I 

To study the anisotropy of the dispersion next to the 
satellite positions we calculate the spin- wave velocities v± 
and W|| perpendicular and parallel to the stripe orienta- 
tion (cf. Fig. ||). For A = 0, where the coupling between 
the domains is switched off, v± is zero and v\\ remains fi- 
nite. With increasing A both velocities increase, v± more 
strongly than Dm . There exists a value A* with isotropic 
velocities, v± = v\\. For p = 4 we find A* sa 0.3. in the 
limit p — > oo both velocities converge to «af as expected, 
for p> lwe find u_l ) ||A'AF ~ 1 ^ 1/P- 111 t ne special case 
A = 1 the velocities are given by 

(24a) 

r_ = — r /' AJ , (24b) 



VAF, 
P 

p-1 



for purely geometric reasons. In this case, all spins are 
interacting in terms of the topology and strength of the 
exchange couplings exactly like in the antiferromagnet. 
The only difference lies in the insertion of strings of holes, 
which effectively stretch the lattice and increase the ve- 
locity by a factor pj {p— 1) in the perpendicular direction. 
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FIG. 3: Spin-wave velocities v±_ and v\\ for vertical stripes 
with spacing p — 4 as a function of A (left panel) and as a 
function of 1/p for different couplings A (right panel; lines are 
a guide to the eye). 

We now focus on the line k = (H, i) containing the 
satellites, along which we plot all p — 1 magnon bands 
in Fig. |] for a variety of p and A. For A < 1 and A > 
1 the bands are separated by gaps. (In this context, 
"gaps" are not necessarily real gaps showing up in the 
density of states, they are apparent gaps along the chosen 
line.) Only for A = 1, the structure seems to consist of 
displaced and intersecting antiferromagnetic bands. The 
value A = 1 is special for the reasons explained above 
which also imply that the band width must coincide with 
the antiferromagnet. The purely geometric effect entails 
just a more complicated band structure. 
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FIG. 4: Band structure for vertical stripes along the (if, 0.5) 
direction with different spacings p and couplings A. Darker 
and larger points correspond to a larger weight of the inelastic 
structure factor. 



To the extent to which our stripe model provides a 
valid description of the magnetic excitations in the ma- 
terials where the n resonance was observed, the resonance 
frequency has to be identified with w(kAF) from the low- 
est magnon band, provided w(kAF) > and the structure 
factor has significant weight. From Fig. ^ one recognizes 
that for A < 1 this is always the acoustical band. On the 
other hand, for A > 1 higher bands may yield a stronger 
resonance (see case p = 3 and A = 2). In Fig. |5| we 
illustrate the dependence of on A and p. 




FIG. 5: The resonance frequency ui-k for vertical stripes as a 
function of A for different spacings p (left) and as a function 
of 1/p for different couplings A (right; lines are a guide to the 
eye). 

For p large enough such that um « i>af and the magnon 
dispersion is roughly linear between the main satellite 



and Icaf, we may estimate 

^tt ~ «af — ■ (25) 
pa 

This estimate becomes exact for small 1/p and represents 
the linear asymptotics in Fig. ^ (right). Deviations grow 
with decreasing p and increasing deviation of A from 1. 

B. Diagonal case 

For diagonal stripes there are more subtle differences 
between even and odd stripe spacings p. Since the basis 
vectors of the magnetic unit cell can be chosen as A' 1 ' 
: 1.1: and A< 2 ) = (p,0) for odd or A( 2 ) = (2p,0) for 
even p (cf. Fig. [j]), we have one hole and p — 1 spins per 
unit cell for odd p and twice the number of holes and 
spins for even p. Like in the vertical case, the number 
of eigenvalues vanishing identically corresponds to the 
number of holes, the number of bands is given by half 
of the number of spins per unit cell, and the bands are 
twofold degenerate. 

All magnetic Bragg peaks are located along the line 
Q = (H,H) with H = j/p for odd and H = j/(2p) for 
even p (cf. Fig. |6|). In the case p = 3 we can calculate 
the dispersion analytically and find 

w(k) = 2 JS{ sin 2 [n(H — K)] + A sin 2 [n(2H + K)\ 

+Asin 2 [7r(!/ + 2X)]} 1/2 . (26) 

Along the k = (H, H) direction this relation simplifies to 

u(H, H) = V2XJS\ sin(37r#)|. (27) 



P = 3 p = 4 p = 5 




H 

FIG. 6: Acoustical band for diagonal stripes with spacings 
p = 3, 4, 5 for A = 1 plotted in analogy to Fig. ^| 

Though the case p = 3 with a single band is the sim- 
plest possible, we find several critical points in the dis- 
persion, which should result in a nontrivial shape of the 
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density of states p(ui) ~ J k <5(w — w(k)). Therefore we cal- 
culate this quantity just to illustrate that even for this 
simplest case p(u>) shows interesting features strongly de- 
pending on the effective coupling A. The numerically cal- 
culated density of states is plotted in Fig. [7] for different 
values of A. The van-Hove singularities are located at 
the energies of the critical points in the dispersion. The 
dependence of these energies on the coupling A is also 
shown in this figure. Due to a finite numerical resolu- 
tion the van-Hove singularities are not resolved if they 
are too close to each other and their precise shape is not 
reproduced, e.g. at the energies of the saddle points p(u>) 
should diverge logarithmically. 
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FIG. 7: Density of states p(u) (height of shaded area in ar- 
bitrary units) for diagonal stripes with spacing p = 3 and 
different couplings A. The thin lines correspond to the ener- 
gies of the critical points in the dispersion. There are up to 
four inequivalent ones, at k = (§,§), k = (|,0), the upper 
band edge with o> ma x and a possible additional critical point. 

Calculating the weight by the structure factor of the 
bands we find the strongest intensity near the zeros of the 
acoustic band at the satellite positions Q = (i±^,|± 

i) for all p. The behavior of the spin- wave velocities v± 
and U|| (cf. Fig. ^) is similar to the vertical case. 

Along the k = (H, H) direction, the acoustical band is 
separated by finite gaps from the optical bands for A 7^ 2. 
For A = 2, the band structure again seems to consist of 
intersecting displaced antiferromagnetic bands. In con- 
trast to the vertical case, the special value of A is now 
2 since for this value the sum of the exchange couplings 
to neighboring spins is as large as in the antiferromag- 
net. However, for diagonal stripes the topology of the 
couplings is different from the antiferromagnet. 

For odd p, the 7r resonance results from the excitation 
of acoustical magnons since the lowest band has a finite 
uj(\i\p) with a relatively strong weight. In contrast, for 
even p the frequency and the weight of the acoustical 
band vanish at kAF- In this case, the tt resonance should 
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FIG. 8: Spin-wave velocities v± and v« for diagonal stripes 
with spacing p = 3 as a function of A (left) and as a function 
of 1/p for different couplings A (right; lines are a guide to the 
eye). 



therefore be ascribed to optical magnons. For A = 2, 
the 7T resonance results from the common edge of the 
acoustical and optical bands (cf. Fig. 0). 
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FIG. 9: Band structure for diagonal stripes along the (H, H) 
direction with different spacings p and couplings A. 



With increasing coupling A the resonance energy in- 
creases. In contrast to vertical stripes, the resonance 
energy remains finite in the limit A — * for even p where 
it arises from an optical band. (cf. Fig. [l0|). Like for the 
vertical case, av decreases with increasing stripe spacing, 
for p> 1 according to uj^ oc I /p. Since resonance comes 
from different bands for even an odd p the tt, is a 
nonmonotonous function of p. For this reason, lOk(p) is 
plotted in Fig. [lC] separately for the two cases. 
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X 1/p 

FIG. 10: Resonance frequency u) m for diagonal stripes as a 
function of A for different spacings p (left) and as a function 
of 1/p for different couplings A (right; lines are a guide to the 
eye distinguishing even and odd p) . 



V. DISCUSSION 

We now discuss our findings in comparison to exper- 
imental data on the spin dynamics, which are obtained 
predominantly from neutron scattering. As a result of 
this comparison we wish to advocate that the simple 
stripe model provides a fair account of the spin dynam- 
ics at not too low energies. At very low energies, spin 
gaps may occur, e.g., due to spin anisotropics (as in non- 
superconducting LNOu and LCOci), due to the coupling 
of spins to the superconducting order parameter (as in su- 
perconducting cuprates, see below), or simply due to the 
absence of antiferromagnetic order (for too small A). Our 
model could straightforwardly be generalized to account 
for the first origin. The inclusion of superconductivity 
would require a major extension. 

In Table [| we have collected basic parameters for var- 
ious undoped compounds setting the fundamental phys- 
ical scales. In Table |j] spin dynamics data for specific 
stripe structures are compiled. 



material 


# layers 


S 


a 


J 


VA¥ 


Refs. 


LNO 


1 


1 


3.8 A 


30 meV 


0.32 eV A 




LCO 


1 


i 

2 


3.8 A 


135 meV 


0.85 eV A 


2e 


27 


YBCO 


2 


1 
2 


3.9 A 


125 meV 


n/a 


i 


c 


BSCCO 


2 


1 

2 


3.8 A 


140 meV 


n/a 


|i 


c 



TABLE I: Basic parameters of the undoped parent com- 
pounds: number of layers in the crystalline unit cell, spin, 
nearest-neighbor spin spacing, nearest-neighbor antiferromag- 
netic exchange coupling, and spin-wave velocity. 



A. LSNO 

We start the comparison with LSNO which displays 
diagonal stripep-,aH|d where integer values of p are par- 
ticularly stableEaO due to a lock-in of the stripes into 
the atomic structure. In this material, static stripes (i.e. 



material 


T c 


P 






Refs. 


LSNO 


OK 


3(d) 


80 meV 


< 28 meV 




16 




LSCO 


« 38K 


4(v) 


n/a 


3.5 meV 


17 


28 


,29 


LSCO 


10K 


6 (v) 


25 meV 


< 1.1 meV 




29 




LSCO 


OK 


« 43 (d) 


7 meV 


meV 




3C 




YBCO 


90K 


5 (v) 


41 meV 


28 meV 


14 


31 


M 


YBCO 


63K 


n/a 


35 meV 


28 meV 




33 




YBCO 


59K 


n/a 


26 meV 


16 meV 




31 




YBCO 


39K 


8(v) 


23 meV 


10 meV 




IS- 




BSCCO 


91K 


n/a 


43 meV 


n/a 


3#5| 


BSCCO 


83K 


n/a 


38 meV 


n/a 




35 





TABLE II: Spin dynamics data for different materials at vari- 
ous doping levels characterized by the critical temperature T c , 
stripe period p and orientation (diagonal/vertical), resonance 
frequency w w , and gap frequency w gap . 



stripes thatj-are visible down to to = Q) are seen at wave 
vectors Q„B3 For p = 3, the spin rhoiamics at higher 
energies has been measured in detail.EEl Similar data are 
also available fprnpninteger p, e.g. p — 3.75.EJ 

Experiments ! 2 3 tr 5 l on undoped material are in agree- 
ment with 2D spin-wave theory for the antiferromagnet 
with J pa 30 meV. This exchange coupling corresponds 
to apJsotropic spin-wave velocity vaf — V&SJa = 0.32 
eVAo since S = 1 and a pa 3. 8 A. This agreement is rea- 
sonably good over a wide energy range u> > 30 meV up to 
the band edge at u> pa 125 meY^at low energies ui < 15 
meV deviations (gaps) appeared due to a uniaxial spin 
anisotropy and weak interlayer couplings. 

The spin dynamics of the stripe.-system was examined 
for p — 3.75 due_tp oxygen dopinged as well as for p = 3 
with Sr doping£3 In the first reduced velocity 

v || ~ 0.6vaf was found in direction parallel to the stripes, 
v± was not resolved. In the second case, the velocity was 
measured in both directions and found to be remarkably 
isotropic and close to the value of the undoped system: 
uy ps 0.30 eVA and v± w 0.35 eVA. The overall shape 
of the magnon dispersion was sinusoidal with an upper 
edge at pa 80 meV. 

In our theory, this sinusoidal shape for p = 3 is well 
reproduced [compare Fig. || and Eq. (p7|)]. The ratio 
u> v /(JS) pa 2.7 is consistent with A pa 0.9. For this value 
of A, vj_ pa vaf and v» pa 0.67«af- Although we find vi\ to 
be smaller than in Ref. the overall agreement is very 
satisfying and provides strong support for our case that 
the spin dynamics can be well understood from a stripe 
model. Small quantitative deviations may be attributed 
to the simplicity of our model using only two types of 
exchange couplings. 

Remarkably, A pa 0.9 implies that the spin exchange 
across a stripe is not much smaller than within an AFM 
domain. It is important to keep in mind that A must 
not be too small to preserve magnetic order. A quantum 
Monte Carlo analysis^ of coupled two-leg ladders (S = 
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1) indicates a quantum phase transition into a disordered 
state at A ~ 0.011. Below this value, stripe order would 
be destroyed by quantum fluctuations. 

Within our approach we can estimate also the two- 
magnon signal accessible by Raman spectroscopy. We 
may compare our single-magnon density p(ui) to the two- 
magnon scattering intensity at frequency 2uj. Certainly, 
this can be made only on a qualitative level, since p was 
calculated neglecting weight factors (which would change 
the shape of spectra but not the frequency of resonances) 
and because linear spin-wave theory does not include in- 
teractions between magnons. Nevertheless, it is instruc- 
tive to compare the outcome from our model for the diaah 
onal case p = 3 with an experiment by Blumberg et aZ.fi3 
on LSNO. In this experiment, two magnetic resonances 
are observed at lu w 4.6J and u> 3J. For A » 0.9 
we expect a singularity in the single-magnon density at 
lu w 2.7JS (see Fig. 0), which would correspond to a 
two-magnon resonance at to 5AJS. If corrections due 
to magnon interactions are modest, the resonance of the 
theory could be idebtified with the upper experimental 
one. Then the resonance at the lower frequency cannot 
be understood. On the other hand, for A not too close 
to 1 the single-band structure for p — 3 would lead to 
several well-separated extrema but contradict the above 
determination of A. In particular, for A < 1, the ad- 
ditional resonance lies above tu^ since it arises from ex- 
trema close to the upper band edge and there is only a 
saddle-point at Icaf- This apparent contradiction might 
be resolved if either interaction corrections are large, ad- 
ditional exchange interaction are important, or the lower 
experimental resonance is of different origin. 



B. Cuprates 

In the present study, we assume the presence of charge 
stripes and evaluate the spin dynamics for a simple 
model. The question of why stripes are formed and how 
stripe formation is related to superconductivity therefore 
cannot be addressed. In particular, the simple spin-only 
model misses the coupling of spin fluctuations to the su- 
perconducting order parameter. Consequently, our anal- 
ysis misses the opening of a spin gap due to supercon- 
ductivity. Therefore, the spin dynamics for to < w gap 
is masked by superconductivity (see Table ||). Never- 
theless, one can expect the stripe-like spin dynamics to 
remain visible in superconducting samples for u > w gap . 

Such a gap has been observed in experiments on LSCO 
(e.g., u igap. 3.5 meV near optimal doping with T c = 
38.5Kjl3IH3'E£l a gap smaller than 1.1 meV for underdoped 
samples with T c = 12K and T c = 25KE2) and on YBCO 
(e.g., Wgap — 10-tneV for a highly underdoped material 



with T r = 39K:la u 



gap 



16 meVJbr a moderately un- 



derdoped material with T c = 59K;I 



"gap 



30 meV for 



near optimal doping with T c = 89Kta). For YBCO there 
is evidencecil for a proportionality between w gap ~ 3.8T C 
which is not far away from the BCS weak-coupling limit 



with Wgap = 3.52T C . 

Furthermore - and more importantly in the present 
context - there is evidence for such a (rough) propor- 
tionality not only between T c and w gap but also. he- , 
tween T c and w w (w, ~ 5T C for underdoped YBCOjf3c3 
~ 5.4T C for under- and overdoped BSCCOE3). From 
our theory, we expect w T to be roughly inversely pro- 
portional to p, see Eq. (|24|) and Fig. [h]. At low doping 
p should be inversely proportional to the doping level 
(x + 28), i.e., u> n should be proportional to the doping 
level. Such a relation was found in a previous theoretical 
study of the Hubbard model,E3 where it was attributed 
to a particle-particle collective mode. Although our ap- 
proach is technically much less involved, it provides an 
alternative explanation which is not in contradiction with 
the previous one since stripe order itself can be consid- 
ered as a collective phenomenon that can be derived from 
the Hubbard model. tlfl 

On the other hand, at larger doping there is no sim- 
ple relation between the doping level and p. In YBCO, 
for example, the charge-transfer mechanism between the 
Cu02 plane and the CuO chains interferes. In LSCO 
it is well .documented that p saturates at p = 4 for 
x > 0.12.E3 Beyond that point (which corresponds to 
optimum dopingo), additional holes may populate the 
antiferromagnetic domains without affecting their period 
p. However, these excess holes may suppress the antifer- 
romagnetic exchange coupling in analogy to holes in the 
spin-glass phase (Ref. 10 reports the corresponding sup- 
pression of the spin- wave velocity). Hence, the effective 
J and, consequently, also may skunk with overdoping 
as seen in experiments on BSCCO.lJ 

For LSCO, so far no direct evidence for a tt resonance 
has been found. This could be simply because the reso- 
nance intensity is expected to be only ~ 10% of the total 



magnetic scattering.Ej However, if the n resonance - in 
the sense of a merger of the incommensurate signals - 
can be attributed to magnons in stripes which are par- 
ticularly well established for LSCO, one definitely should 
expect such a resonance. For underdoped LSCO (p = 6, 
T c = 25K) there is evidence for = 25 meV (where 
incommensurate response becomes commensurate) BB A 
similar signal was observed at even lower ^doping in the 
spin-glass phase (w T = 7 meV for p w 43) Ei 

Like for LSNO, we may use the values of J, p, and u)„ 
to estimate A for the cuprates. Eop-YBCO with J = 125 
meV, p = 5, and = 41 meV,EJo we obtain A ~ 0.07 
from the left panel of Fig. |5| If we take—/ = 135 meV and 
u) n — 25 meV for LSCO with p = 6,E£l an even smaller 
value A ~ 0.04 is found. 

^From this result we may predict where the resonance 
should be expected in LSCO near optimal doping (p = 
4). For J = 135 meV and A = 0.04 - 0.07 we find tu^ » 
40 — 52 meV. While the resulting values for have a 
certain spread, they suggest that the resonance frequency 
should be at least as large as in optimally doped YBCO. 

In the experiments known to us, the considered energy 
range was simply too small to detect the resonance for 
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optimallydoped LSCO: lo < 6 meV mj§[ w < 10 meV 
in Refs. (45|]4q , u> < 16 meV in Ref. |17[ However, from 
pulsed neutron scattering evidence has been found for 
a broad peak in the moawjBturn-integrated susceptibility 
between 40 and 70 meVO'Ell which could be ascribed to 
the 7r resonance. 

Apparently, A seems to be significantly smaller in the 
cuprates than in the nickelates. At the same time, S is 
smaller (although J is larger). Therefore one may wonder 
whether static magnetic stripe order is already destroyed 
by quantum fluctuations without invoking competing or- 
ders leading to a gap. For S = \ the coupling.meeds to 
satisfy A > 0.3 to stabilize spin order-Jor p = 3,EiE3 while 
for p = 4 a finite A > is sufficient For p = 5 (as for 
every odd p) one again expects a finite critical A. If the 
intersptripe coupling is below this value, the presence of 
a spin gap can be understood also within the spin-only 
model. 



C. Conclusion 

In summary, we find that the spin fluctuations of 
stripes can provide a simple and valuable description of 
the dynamics observed in high-T c compounds and related 
materials. Already our minimalistic spin-only model pro- 
vides an accurate account of experiments on LSNO and 
possibly also a unifying framework for incommensurate 
response and the it resonance in the cupratea. While 
such a framework has been suggested recentlytJ it is an- 
alyzed and evaluated here for the most transparent case 
of integer periods p. Our results unravel the evolution of 
the band structure with p for diagonal and vertical stripe 
configurations. Likewise, we have explicitly determined 
the dependence of characteristic spin- wave velocities and 
of the resonance frequency on p and A. Thereby, we pos- 
tulate that the n resonance reflects the magnon frequency 
u>n of the lowest lying band with nonvanishing weight. In 
particular, uj^ was found to be roughly inversely propor- 



tional to p in agreement with experiments. 

Hopefully, future experiments can provide more direct 
evidence for the it resonance also in LSCO. This would 
also relax the controversial question, whether spin exci- 
tations in LSCO and YBCO are analogous! 7 ] or not a If 
stripe magnons indeed explain the spin dynamics at in- 
termediate energies, as we expect, they would provide a 
unifying framework for understanding the spin dynamics 
above the gap scale. Then the stripe physics would be 
also of great importance as basement for superconduc- 
tivity as low-energy phenomenon. 

Naturally, several aspects remain unexplained by our 
minimalistic theory. For example, our model cannot be 
expected to explain why the magnetic incommensurabil- 
ity disappears at XL in YBCOEfl while charge order is 
visible up to 300K.E3 Probably this is a question to the 
stripe-forming mechanism and to a possible coupling be- 
tween the order parameters for stripe order and super- 
conductivity. In LSCO, the vicinity of soft phonons and 
structural instabilities may help to stabilize stripes at 
temperatures above the superconducting transition. 

For future studies it would be interesting to include 
effects of the bilayer coupling present in YBCO and 
BSCCO, of the weak 3D coupling present in all mate- 
rials, as well as spin anisotropy, more complicated spin 
interactions (e.g. four-spin cyclic exchangeEJcj) , excita- 
tions beyond spin waves (e.g. double-spin excitationsej) , 
mobility of spins, and effects of disorder, to name just a 
few. 
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